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The results of a relaxation-mode analysis are presented for two cases of trapped-defect 
relaxation in the NaCl structure, in which both defects occupy the same type of site (e.g., 
impurity divalent ion and trapped vacancy), or in which they occupy the two different 
types of sites (vacancy pair). The relaxation analysis is presented in the form of a set of 
basis vectors in occupation-probability space and a set of secular equations. Solutions to the 
equations provide the relaxation rates and also the coefficients in the linear combinations of 
the basis vectors which constitute the relaxation modes. 

Calculations of the relaxation rates and contributions to the polarizability of the various 
modes for a three-shell model with jump frequencies chosen to represent the relaxation of 
an impurity-ion vacancy pair in NaCl(Mn) have confirmed the results of Lozovskii. Even 
though the third shell makes a significant and even large; contribution to the process, one 
mode dominates and the relaxation as seen in a-C measurements would take place with es- 
sentially a single relaxation time. The d-c techniques of Dreyfus are sensitive enough to 
detect more than one of these; relaxations at low temperatures, but the slowest of these will 
always be the dominant one. Any relaxation slower than the major one must be ascribed to 
some other mechanism. 

1. Introduction 

Point defects in crystals give rise to mechanical and elastic relaxation in a variety of ways. 
In particular, because of their interaction with each other, pairs of them can become mutually 
trapped, forming centers which can contribute to the response to applied mechanical or electrical 
stresses. If one defect is considered as the trap, then the trapped defect can occupy various 
lattice sites in the neighborhood of the trap. These accessible sites may be restricted to those 
in the nearest neighbor shell to the trap because of very strong interaction, or they may include 
sites in shells further away from the trap. Under the influence of the applied stress, some of 
the accessible sites become preferred over others that are equivalent in the stress free crystal, 
and the redistribution of the trapped defects according to the new scheme constitutes the 
relaxation process. 

The rate equations governing this redistribution process constitute an eigenvalue problem, 
for which group representation theory combined with a knowledge of the symmetry around the 
trap provide a powerful tool for solution [1-3]. 1 The results of this treatment are expressed in 
terms of relaxation modes, which represent certain distributions of the trapped defects among 
the accessible sites. Each of these distributions relaxes as a unit with the same exponential time 
function, exhibiting a single relaxation time, governing the occupation probability of all the 
sites. 

Haven and van San ten [1, 2] applied the group theory technique to the case of cation 
vacancies trapped near divalent cation impurities in the NaCl structure, with only the nearest 
neighbor sites included. Wachtman [3] showed that the relaxation modes of the undriven 
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situation were also solutions for the driven case in the limit of small applied stresses, and 
applied the technique to the relaxation of oxygen vacancies trapped at divalent cation im- 
purities (Ca ions) in Th0 2 , again restricting the vacancies to nearest neighbor sites. Bhaga- 
vantam and Pantulu [4] have summarized this approach for the case of nearest neighbor sites 
and have emphasized the similarity of the symmetry arguments to those arising in the treat- 
ment of infrared and Raman spectra. Franklin [5] showed that the technique could be ex- 
tended to the case in which additional shells, further out than the nearest-neighbor shell, were 
accessible, and gave solutions for the case of cation vacancies trapped at divalent cation 
impurities in the NaCl structure, including two shells of accessible sites. 

In considering the dielectric and mechanical behavior of these centers, the more distant 
shells of sites must, at least in principle, be included. It is expected that, at least beyond the 
first two shells, the binding energy of the sites will decrease with distance from the center, 
and the probability of occupation will fall off in the more distant sites. However, there are 
more sites available in the more distant shells. In the dielectric case, also, the more distant 
sites will make a larger contribution through their larger dipole moments, so that even a 
relatively small occupation probability could be important. The influence of the more dis- 
tant shells is expected to be important when the differences in binding energies of the shells 
are not too large, and when the temperatures are high. 

The experimental data are best defined in the case of dielectric relaxation in the alkali 
halides containing divalent impurities [6-10]. When care is taken to avoid precipitation of 
the divalent impurity, the frequency spectrum of the relaxation ascribable to these centers 
can be described as a typical Debye peak, with but a single relaxation time. Lidiard [11] has 
shown that this is to be expected if the first two shells of neighbors are included in the calcu- 
lation, and reasonable values assigned to the probabilities for jumps of the vacancy between 
and within shells. Recently, Dreyfus [12], using a d-c technique of considerably greater time 
resolution than found in the usual a-c methods, was able to separate two distinct relaxations 
for which he could account using Lidiard 's treatment, and another, slower, relaxation which 
he has tentatively ascribed to the third and more distant shells. 

Lozovskii [13] has made a relaxation mode analysis for this dielectric relaxation including 
three shells of sites, but considering two rather special cases. For binding energies of the 
various shells, he chose (1) purely Coulombic values arising from the electrostatic interaction 
between point charges in a dielectric medium, and (2) the same value for all shells. He cal- 
culated both the relaxation times relative to the jump time in the first shell and the relative 
contributions to the polariz ability from each mode. His results indicate that there should be 
four active modes, but that one will predominate, so that a close approximation to a Debye 
peak should result even when the third shell is included. Furthermore, the dominant relaxa- 
tion in both cases was the slowest one. Dreyfus "fast" and "very fast" polarizations can be 
ascribed to the dominant mode and a superposition of the other three; but since the dominant 
mode should be the slowest, Lozovskii's result suggests that Dreyfus' "slow" polarization 
should not be ascribed to the outer shells. 

In part 1 of the present paper, we extend the mode analysis to the third shell of neighbors 
for two cases in the NaCl structure. The Matched Point Defect Pair, in which the trap and 
the moving defect are both on lattice sites of one kind, as for a cation (anion) vacancy trapped 
near a cation (anion) impurity; and the Unmatched Point Defect Pair, in which the trap is 
on one kind of site and the moving defect is on the other, as in a vacancy pair, where either 
vacancy can be considered as the trapped and the other as the moving defect. This analysis 
includes all possible relaxation modes under all possible stresses, and includes a discussion of 
the determination of which modes are excited by a given applied stress. In part 2, we present 
calculations similar to those of Lozovskii for the relaxation times and polarizabilities for the 
dielectric relaxation of the three-shell model of a Matched Pair. These calculations have 
been made with a set of binding energies chosen to provide as realistic as possible a description 
of the situation in NaCl (Mn), and the results are compared to the available data. 
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2. Relaxation Modes of Pairs of Point Defects in the NaCl Structure 
2.1. Matched Defect Pair (Impurity Ion — Vacancy) 

This calculation is a straightforward extension of that reported previously [5] for the two- 
shell model. The modes for the three-shell model will be linear combinations of the previous 
ones together with new contributions from the third shell. Because of the complexity of the 
situation, with 42 sites available to the moving defect, it is preferable to report the modes in 
tabular form rather than graphically as done before. A 42-fold vector space is defined by the 
occupation probabilities of the accessible sites. Any distribution is then expressed as a vector 
p in this space. The rate equations are expressible in matrix form, involving a matrix C whose 
elements are c tj , the negative probabilities per unit time that the moving defect will jump from 
the ith site to the jth site, and cu, the sum of all jump probabilities out of the ith site. The 
solutions to the rate equations are then expressible as sums of the eigenvectors of the matrix 
C, 

P=\ ^i + S Mt \Y) u n) , where /i=/ i0 exp [—\ t t]. 

In this equation, the u f are the normalized eigenvectors of C, U\ is the eigenvector for which the 
eigenvalue (jump rate) is zero and which represents the equilibrium distribution, and the 
U\ k are the components of Si. A^is the total number of accessible sites, the X* are the eigenvalues 
(jump rates, or reciprocal relaxation times), and the j i{) are determined by the initial distribution. 

In the irreducible representation technique, the eigenvectors % % are expressed in terms of 
a set of normalized basis vectors v V 'To each irreducible representation of the point group of 
the relaxing center in the vector space of site occupation probabilities corresponds a subspace, 
spanned by one set of these basis vectors, and from which an equal number of eigenvectors can 
be constructed. If a transformation matrix is constructed of all of these vectors as columns, it 
can then he used to transform the matrix C to reduced form. The problem has then been re- 
duced to a set of secular equations, one for each irreducible representation whose solutions pro- 
vide the eigenvalues and the coefficients in the linear combinations of the rectors v ^ that form 
the eigenvectors u t . 

In table 2.1 below, we list the 42 basis vectors used for this problem. The first column 
numbers the 42 sites in the first three shells of neighbors, and the second column identifies their 
lattice position with respect to the impurity atom, relative to the cation-anion separation. In 
the succeeding columns are listed the vectors used as basis vectors, arranged by irreducible rep- 
resentations, which are identified at the top of the columns. As it stands, this table is also the 
transformation matrix, except that for convenience of presentation, the basis vectors have not 
been normalized. The irreducible representations are identified in the notation of Eyring, 
Walter, and Kimball [14]. 

The secular equations corresponding to each irreducible representation are presented next. 
It is understood that in solving these secular equations, the normalized basis vectors must be 
used. The w tj are the probabilities per unit time that the vacancy in a site in the ith shell will 
jump into a neighboring one in the /th shell. Only jumps from a site to one of its 12 nearest 
neighbors are allowed. For the interchange of the impurity ion and a vacancy on a first shell 
site, w f is used. Analytical solutions for the eigenvalues and for the coefficients in the linear 
combinations of the basis vectors to form the eigenvectors have not been obtained for the E g> 
T iu and the T 2g representations. 
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2.2. Unmatched Point Defect Pair (Vacancy Pair) 

The method of calculation is similar to that for the matched point defect pair. Each 
configuration can be changed into another accessible one by jump of either the cation or the 
anion defect. Hence, the jump probabilities always occur as the sum of that for the cation 
and that for the anion defect, and in the notation used here, only these sums will be used. It 
should also be noted that there are two different types of jumps from a site in the third shell to 
another in the same shell, distinguished by the path taken. For example, the sites in the third 
shell are designated by (2,1,0) and its permutations, including negative signs. A jump can 
only be made to another site in this shell for which the x,y,z differ by ± 1. If the zero stays 
fixed, there is only one such jump, and the defect remains in a (100) plane as it moves. During 
the jump, it is screened from the other defect by an ion of opposite sign, in the same (100) 
plane. The other two possible jumps are with the two remaining fixed, for which the path is 
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not screened. Presumably, the activation energies for these two jumps will not be the same, 
and so they will have different jump probabilities. These will be designated by w^ for the 
screened and W33 for the unscreened jump. With these changes, the notation as used for the 
previous example will be used to designate the jump probabilities. 

The basis vectors are given in table 2.2, again using the cation-anion separation lattice 
parameter as the repeat distance in designating the sites. The secular equations are given 
below from which the eigenvalues and coefficients in the linear combinations of basis vectors 
constituting the eigenvectors can be derived as before. 
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2.3. Excitation of the Modes by Applied Stresses 

While the analysis above reveals the modes by which a relaxation process can take place it 
does not provide any way of deciding which modes will be excited by a given applied stress. 
This question can, however, be answered in the following way. Consider the system under a 
constant stress of the specified sort. The influence of the stress will be to lower the symmetry 
of the crystal to that common to the applied stress and the undisturbed crystal (Curie's law [15]). 
Thus, from a knowledge of the symmetry of the applied stress and the undisturbed point group 
of the defect center under consideration, the point group under stress can be determined. It 
will be a subgroup of the unstressed point group. In the presence of the stress, the system will 
adopt the equilibrium configuration appropriate to the new point group, and represented by the 
new A lg mode, for which the relaxation rate is zero. When the stress is removed, the system will 
begin to relax, but now with the modes appropriate to the undisturbed point group. Hence, 
the problem of deciding which modes will be excited is the problem of expressing the equilibrium 
mode under the stressed point group in terms of the unstressed modes. 

The irreducible representations of the unstressed group each correspond to a subspace of 
the total vector space, such that these subspaces are mutually orthogonal. Under the stressed 
group, the equilibrium mode corresponds to a single vector, which must of necessity be a linear 
combination of mutually orthogonal components drawn from the appropriate subspaces 
identified with the irreducible representations of the unstressed group. Since each of these 
components, in turn, can be expressed as a linear combination of the basis vectors of the 
corresponding unstressed irreducible representation, then the stressed equilibrium mode must 
be expressible as a linear combination of the unstressed modes whose corresponding irreducible 
representations have as a component the stressed irreducible representation A ig . Hence, the 
problem is to examine the irreducible representations of the unstressed group and determine 
which of these has as a component the stressed A lg irreducible representation. 

This determination is readily done with the character table. Since the irreducible repre- 
sentations of a group are also representations, although not necessarily irreducible, of its sub- 
groups, the character table of the irreducible representations of the full group can be consulted, 
and the examination restricted to the entries under the elements of the subgroup. For a given 
irreducible representation of the full group, these constitute the characters of that representation 
of the subgroup, and by inspection and comparison to the character table for the subgroup, the 
components among the irreducible representations of the subgroup determined. In this way, 
all of the irreducible representations of the full group can be found that have as a component 
the Aig representation of the subgroup. The corresponding relaxation modes of the unstressed 
system are those that will be excited by the stress that produced the subgroup. 

As an example, consider the case of the matched point defect pair in NaCl, and let the 
stress be an electric field applied along one cube axis. The unstressed group is the full cubic 
group, m3m in the international notation; while under the field, this is reduced to 4mm. In- 
spection of the character tables reveals that the unstressed irreducible representations which 
have the A x irreducible representation of the stressed subgroup as a component are the A Xg , the 
E gy and the T lu . On the other hand, a simple tension or compression along a cube axis reduces 
the point group to 4/mmm, and results in excitation of the A ig and the E g . This result can be 
interpreted as follows. The A lg mode will always be excited in principle, because a stress will 
never shift the distribution completely away from the equilibrium distribution, and there 
must be a term to describe what is left. The E g mode represents the departure from equilibrium 
that provides the symmetric extension or contraction along the cube axis produced by the 
tension or compression. Under the axial electric field, this mode is also excited and corresponds 
to the electrostriction. 

There is a rather simple way of determining the number of independent relaxation times 
to expect for the excited modes that parallels the above development. Under the applied stress, 
in the new point group, the sites will have been further subdivided into independent and non- 
interacting subspaces, analogous to the various shells of neighbors of the unstressed crystal. 
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The Aig mode in the stressed point group will be the sum of vectors, one from each subspace, 
each formed simply by assigning unity to each site in the subspace. These vectors are the 
mutually orthogonal components making up the new equilibrium as discussed above. The 
number of distinct relaxation times observable in principle under the applied stress will be just 
the number of these independent subspaces produced with the proviso that one of these relaxa- 
tion times will be infinite, corresponding to the equilibrium mode, and the others finite. 

This treatment can, of course, be generalized. A listing of all subgroups reached by the 
application of any homogeneous second rank tensor stress to any crystal was given by Wacht- 
man and Peiser [16]. In table 2.3, we complete this tabulation by listing all subgroups reached 
by the application of a vector stress, such as an electric field. Koptsik [17] has tabulated the 
subgroups reached by various types of influence, but his tables apparently do not apply to the 
present case. The subgroups which can be reached by application of an electric field have 
previously been listed by Vlokh and Zheludev, [18] but the meaning of their table for the cubic 
groups is hard to interpret, at least in the English translation. Table 2.3 also gives some 
additional information. The subgroups are classified as to whether they are normal subgroups 
(sometimes termed self -con jugate or invariant subgroups) of the starting point group. This 
classification shows that the cryptosymmetry theory of Niggli [19] and Wittke [20] apparently 
cannot be applied directly to the treatment of relaxation modes. In Niggli's theory, each atom 
is considered to have some additional property, such as magnetic moment, and the symmetry 
group possessed by the crystal when this additional property is taken into account must be a 
subgroup of the group obtained when this additional property is ignored. One might attempt 
to take the probability of occupation by a vacancy as the additional property and discuss 
relaxation modes in terms of cryptosymmetry. However, in Xiggli's theory, only normal 
subgroups are permitted, and table 2.3 shows that many of the subgroups which can be reached 
by the action of an electric field are not normal subgroups. It thus appears that some modi- 
fication or extension of the concept of cryptosymmetry would be required if it were applied 
to relaxation modes. 
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Table 2.3 also shows that some groups can occur as their own subgroups under the action 
of electric field wmile others cannot. This distinguishes between cases in which it is possible 
to apply a suitably oriented electric field which will excite no relaxation mode (except, of 
course, the trivial mode) and other cases in which any field no matter what its orientation will 
excite a nontrivial mode. It should be noted that if a purely arbitrary direction of the stress 
axis be chosen, it can always be arranged that no symmetry elements exist in common with 
the unstressed point group and the stress, so that the symmetry is reduced to the identity 
operation 1, and all modes are excited. The application of a homogeneous tensor stress to a 
centrosymmetric crystal can reach T, so that any centrosymmetric mode can be excited. These 
results have their parallel in the expansion of the thermodynamic potential of a crystal in 
terms of elastic, dielectric, piezoelectric, electrostrictive, etc., coefficients, including higher order 
terms. 

The ability of a suitably oriented electric field to excite all relaxation modes does not 
conflict with the classification of modes by Haven and van Santen [1, 2] into electrically ac- 
tive, mechanically active, etc., provided that their classification is taken to mean the modes 
which are most strongly excited by the appropriate influence with the understanding that 
other modes may be less strongly excited. One can understand how this occurs by consider- 
ing the effect of an electric field upon an electrically active mode and upon a mechanically 
active mode. The effect of a static electric field is to shift the energy levels of the sites. These 
energy shifts, E h affect the occupation probabilities exponentially through their appearance 
in Boltzmann factors; each exponential can be expanded in a power series and only the first 
few terms kept if the energy shift is small compared to JcT. The zeroth order terms are each 
multiplied by the shift in occupation probability for the corresponding site and then summed 
to give the change in potential energy. For an electrically active mode (one with a dipole 
moment) the sum will not be zero, and the mode will be excited to first order in the field; for 
a mechanically active mode, the sum must be exactly zero because the mode is centrosymmetric 
so that the mode will not be excited to first order in the field. However, the second order 
terms from the power series will all be of the same sign and when multiplied by a centrosym- 
metric pattern of occupation probability and summed, these will give a nonzero result so that 
a mechanical mode can be excited in second order by an electric field. This line of argument 
leads to the following rule for the excitation of relaxation modes in terms of their lowest order 
nonzero moments : Modes having a dipole moment can be excited in the first order by electric 
field ; those with a quadrupole moment can be excited in the first order by mechanical stress 
or in the second order by electric field; those having an octapole moment can be excited in 
the third order by electric field; those having a 16-pole moment can be excited in the second 
order by mechanical stress and in the fourth order by electric field. 

Chang [22] has recently tabulated the mechanical modes (in the above sense) for paired 
point defects in a number of structures with m3m point group symmetry. 



3. Dielectric Relaxation in NaCl Containing Divalent Impurity Ions 

In order to assess as realistically as possible the influence of the third shell on the dielectric 
behavior of a real system, we have calculated the dielectric relaxation times and polarizabilities 
using jump frequencies estimated to correspond to a vacancy trapped near a Mn ++ ion in 
NaCl, and employing the results of the present paper. The jump frequencies were chosen 
as follows. It was assumed that the jump frequency from the ith to the jth site could be 
given by w n exp [(Ej—E^fekT], where the Ej are the binding energies of the various sites, 
and the w u is the jump frequency in the first shell. Furthermore, following Haven [21], w 33 
was taken to be % w u and Wi was neglected. For the binding energies, the value of 0.41 eV 
given by Watkins [23] for the first shell was used, and the value of 0.38, derived from his deter- 
mination of 0.034 eV as the difference in binding energies between the first and second shells, 
was chosen for the latter. Coulombic energies were plotted for the first nine shells, and a 
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Figure 1. Relative polarizabilities for the matched point 
defect pair, in units of (ea) 2 /kT where e is the charge on 
each defect and a is the nearest neighbor distance, for 
models allowing from 1 to 9 shells of neighbors. 
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smooth curve drawn through Wat kins two values and joining the Coulombic curve at and 
above the fifth shell. Temperatures of —46° and 227° C were chosen to match the experimental 
ranges of Dreyfus [12] and Watkins [23] respectively. 

With these values for the binding energies, an estimation of the polarizability per center 
could be made for models allowing access to various numbers of shells of neighbors. The 
relative polarizabilities per center, in units of (ea) 2 /kT where e is the charge per defect and 
a is the distance between nearest neighors, are shown in figure 1 for the two temperatures. 
At —40 °C, little error is made in stopping at the third shell, but at 227 °C, the magnitude 
of the polarizability is at least three times that estimated From a single shell of neighbors. 
The influence of outer shells certainly cannot be neglected, at least at the higher temperatures. 

The contribution of each of the electrically active 4 modes to the total polarizability is given 
in table 3.1 . The outstanding feature is the dominance of the slowest mode. This contributes 
96 percent of the total polarizability at both —46 °C and 227 °C, and leads to the conclusion 
that a Debye peak would result from an a-c study of the dielectric relaxation with an apparent 
single relaxation time, even if the influence of the third shell were included. 

In the present calculation, the third shell contributed 46 percent of the polarizability 
at 227 °C, and 10 percent at —46 °C. It appears that the dominance of one mode in the 
polarizability may be a general feature of these models, regardless of the number of shells 
included. 

The effect of temperature upon the relaxation can be seen from the figures in table 3.1. 

Table 3.1. Dielectric relaxation of S-shell model 
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4. Summary 

The results of a relaxation mode analysis are presented for two cases of trapped-defect 
i elaxation in the NaCl structure, in which both defects occupy the same type of site (e.g., 
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impurity divalent ion and trapped vacancy), or in which they occupy the two different types 
of sites (vacancy pair). The relaxation analysis is presented in the form of a set of basis 
vectors in occupation-probability space and a set of secular equations. Solutions to the equa- 
tions provide the relaxation rates and also the coefficients in the linear combinations of the 
basis vectors which constitute the relaxation modes. 

The modes excited by the application of a given stress can be deduced from the subgroup 
of the full symmetry group of the crystal reached by action of the stress. Under this sub- 
group, some of the group representations that were irreducible under the full group become 
reducible. Some of them contain the A lg irreducible representation of the subgroup as a com- 
ponent. The corresponding modes are then the ones that are excited by the stress. 

Calculations of the relaxation rates and contributions to the polarizability of the various 
modes for a three-shell model with jump frequencies chosen to represent the relaxation of an 
impurity-ion vacancy pair in NaCl (Mn) have confirmed the results of Lozovskii [13]. Even 
though the third shell makes a significant and even large contribution to the process, one mode 
dominates and the relaxation as seen in a-c measurements would take place with essentially a 
single relaxation time. The d-c techniques of Dreyfus [12] are sensitive enough to detect more 
than one of these relaxations at low temperatures, but the slowest of these will always be the 
dominant one. Any relaxation slower than the major one must be ascribed to some other 
mechanism. 



The authors thank H. S. Peiser for help with table 2.3. 
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